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We evaluate the percolation threshold values for a realistic model of continuum segregated systems, 
where random spherical inclusions forbid the percolating objects, modellized by hard-core spherical 
particles surrounded by penetrable shells, to occupy large regions inside the composite. We find 
that the percolation threshold is generally a non-monotonous function of segregation, and that an 
optimal (i. e., minimum) critical concentration exists well before maximum segregation is reached. 
We interpret this feature as originating from a competition between reduced available volume effects 
and enhanced concentrations needed to ensure percolation in the highly segregated regime. The 
relevance with existing segregated materials is discussed. 
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The percolation threshold of a two-phase heteroge- 
neous system denotes the critical concentration at which 
global (long-range) connectivity of one phase is first es- 
tablished, and is accompanied by a sudden transition 
of the effective properties of the whole system, such as 
the conductivity in conductor-insulator composites or the 
permeability in viscous fluids flowing through porous me- 
dia [H, m . Unlike the universal (or quasi- universal) be- 
havior of the critical exponents characterizing the per- 
colative transition, the value of the percolation threshold 
is a function of several variables such as the shape of the 
percolating objects, their orientation and size dispersion, 
their possible interactions and the microstructure in gen- 
eral [3|. 

Of fundamental importance for several technological 
applications is the possibility of exploiting such a multi- 
variable dependence to tailor the percolation threshold. 
In particular, an issue of great interest concerns the prob- 
lem of lowering the percolation threshold, so to have long 
range connectivity of the percolating phase with the min- 
imum possible critical concentration. This is the case 
when, for objects dispersed in a continuous medium, one 
wishes to exploit the properties of the percolating ele- 
ments, but still preserving those of the host medium. For 
example, low conducting filler amounts in a conductor- 
insulator composite permit to obtain an adequate level 
of electrical conductivity with mechanical properties of 
the composite being basically unaltered with respect to 
those of the pristine insulating phase. 

There exist two main strategies to lower the perco- 
lation threshold through the manipulation of the mi- 
crostructure of heterogeneous composites. In one of such 
methods, one exploits the large excluded volume Vgx 
of particle fillers with large aspect-ratios, such as rods 
and/or disks dispersed in a three dimensional continuum 
medium [4], whose critical concentrations, being propor- 
tional to l/^ex, can be made extremely small for suffi- 
ciently large aspect-ratios [5]. Such percolation thresh- 
old lowering has been studied in detail for several parti- 
cle shapes and inter-particle interactions, and is now well 
documented ml. 



A second strategy to lower the percolation threshold is 
obtained by forbidding the percolating objects to occupy 
large (compared to the particle size) volumes inside the 
material, so to give rise to a segregated structure like the 
one shown in Fig.[T](b). In practice, this can be achieved 
when elements of two (mutually impenetrable) species 
have different sizes and percolation is established by the 
smaller elements. Typical examples of segregated sys- 
tems are conductor-insulator composites where the size 
of the conducting particles is much smaller than that of 
the insulating regions 0, S, 0] , which display critical con- 
centrations of a few percent or lower. 

Despite that such low percolation thresholds are qual- 
itatively understood by the reduced available volume for 
arranging the conducting particles, very few studies exist 
on segregated percolation in the continuum [10], while 
the vast majority of studies are limited to lattice rep- 
resentations of the segregated structure [1, [ll|, which 
provide only a partial understanding of the percolation 
properties of segregated systems. 

In this Letter we consider a realistic continuum model 
of segregated percolation, primarily aimed at describ- 
ing the microstructure of segregated conductor-insulator 
composites, but general enough to represent also other 
structurally similar systems such as particle-laden foams 
[l2| or filled asphaltene matrices [13]. We show that, 
by varying the degree of segregation of the system, the 
percolation threshold is generally not a monotonous de- 
creasing function of segregation, as suggested by earlier 
studies [1, [13, [Hi, but rather it displays a minimum be- 



fore maximum segregation is reached. Hence, the optimal 
percolation threshold does not necessarily coincide with 
the most segregated structure, leading to a more complex 
phenomenology than previously thought. 

We model a continuum segregated composite by con- 
sidering one kind of impenetrable spherical particles of 
diameter <Ji, which may refer to the conducting objects 
in a conductor-insulator composite, and a second kind of 
(insulating) spherical particles with diameter (72 > ai, 
which we allow to penetrate each other. Furthermore, to 
generate segregation, we assume that the two species of 
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FIG. 1: (Color) Percolating cluster of the conducting phase 
for (a) the homogeneous case 02 = and for (b) the segre- 
gated regime with G^jcri = 12 and 02 = 0.89. The conducting 
particles are plotted together with their penetrable shells with 
d = ai. The color map defines the values of the connectivity 
number k for each particle (see text). 



particles are mutually impenetrable, and that the voids 
left over from the two kind of particles are filled by the 
second (i.e., insulating) phase. Finally, the connectivity 
criterion for the conducting phase is defined by intro- 
ducing a penetrable shell of thickness d/2 surrounding 
each conducting sphere, so that two given particles are 
connected if their penetrable shells overlap. This model 
represents a rather faithful description of real segregated 
composites, such as the Ru02-glass systems [1, where 
thermal treatments on mixtures of Ru02 and glassy 
grains lead to composites made of conducting Ru02 par- 
ticles dispersed in a continuum insulating glassy phase. 
Segregation is induced by the larger size of the original 
glassy grains compared to that of the conducting parti- 
cles. Furthermore, in this and other similar classes of 



composites, electrical transport is given by direct tunnel- 
ing or hopping processes, defining a characteristic length, 
represented by d in our model, below which two con- 
ducting particles are electrically connected. Finally, it 
is worth to mention that the model introduced here is 
relevant also for studying transport of macromolecules in 
disordered porous media [1^, where d represents in this 
case the size of a test macromolecule and (72 / cri the pore 
size ratio in a bi-dispersed porous medium. 

In our numerical simulations, the system described 
above is generated by first placing randomly the insu- 
lating spheres in a cube of edge length L with a given 
number density p2 = N2/L^, where N2 is the num- 
ber of spheres. The corresponding volume fraction for 
L ^ 00 is 02 = 1 — ex.-p{—V2p2)^ where V2 = Traf/O is 
the volume of a single insulating sphere. In a second 
step, Ni conducting (and impenetrable) particles of di- 
ameter (Ji and number density pi = Ni/L^ are added 
in the remaining void space and a Metropolis algorithm 
is used to attain equilibrium In the following, for 
the conducting phase, we use the reduced concentration 
variable r]i = pi7r/6(cri + d)^. Examples of the result- 
ing spatial distributions of the conducting particles are 
shown in Fig. W^a) for the homogeneous case (02 = 0) 
and in Fig. [Tfb) for a segregated one {(J2/(Ti = 12 with 
02 = 0.89) for 7^1 values close to their respective percola- 
tion thresholds (see below). 

To obtain the critical density rj^ for infinitely large sys- 
tem sizes we follow a standard finite-size scaling method. 
Namely, for given values of (12/(71, and p2, and by us- 
ing a modified Hoshen-Kopelman algorithm [15], we cal- 
culate as a function of rji and L the probability P{r]i^ L) 
that a cluster of phase 1 spans the system in a given direc- 
tion, with periodic boundary conditions in the other two 
directions. The critical density r]i{L) for finite L is then 
extracted from the condition P{r]i^L) = 1/2 [l^], and 
r]i follows from the scaling relation r]i{L) — rji oc 
where u is the correlation length exponent obtained from 
the width of the transition. We considered 8 different 
system sizes ranging from L = 16 with Ng = 1500 re- 
alizations to L = 60 {Ng = 100) for (J2/(Ti = 1 and 
from L = 60 {Ns = 200) to L = 140 {Ns = 100) for 
cr2/cri = 12. 20 values of r]i were typically used to fit 
P(?7i,L) with an appropriate function. In this way, for 
most of the cases studied, the calculated v values were 
well within 5% of the universal value 0.88 [l|. 

Typical spanning probability results are reported in 
Fig.O where we plot P{r]i,L) for (12/(71 = 4, d = (Ji, and 
for two values of 02 with few different system sizes. As 
it is clear from the figure, compared to the homogeneous 
case 02 = 0, the spanning probability transition for 02 7^ 
gets shifted to lower values of 771, indicating that the 
percolation threshold is reduced by segregation. This is 
confirmed by the scaling analysis described above, which 
gives r]i = 0.3203 ± 0.0003 for 02 = 0, which is in very 
good accord with Ref.[i3,[i3, and r]^ = 0.1821 ± 0.0004 
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FIG. 2: (Color online) Spanning probability as a function 
of r]i for few values of the system linear size L and for two 
different values of insulating phase volume fraction 02. The 
penetrability length is d = ai and a2/(Ji = 4. 
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FIG. 3: (Color online) Percolation threshold values r/i as a 
function of the volume fraction 02 of the insulating spheres for 
(a) d = ai and several value of G^jcJi and (b) G^jcFi =4 and 
few values of d. The vertical dashed lines are lower bounds 
of the maximum segregation obtained from Eq.([T]), while the 
dotted lines are from Eq.(|2]). 

for 02 = 0.65. 

Although the reduction of r]\ shown in Fig. [2] has to 
be expected on the basis of reduced available volume ar- 
guments, we find that, actually, this is not always the 
case, and that r]\ is generally a non- monotonous function 
of 02 for fixed (12/ (Ji. This is shown in Fig. [3](a) where 
rii is plotted as a function of 02 for d = cfi and for sev- 
eral values of (T2I (y\^ and in Fig. [3l[b) where (72 /ai = 4 
and d is varied. For all cases studied, on enhancing 02 
from the homogeneous case 02 = 0, the behavior of the 
percolation threshold is characterized by an initial linear 
decrease of r^f, followed by a minimum at a particular 
value of 02 which depends upon G2 / cri and and a final 
increase well before maximum segregation is reached at 
02- Lower bounds of 0^ are plotted in Fig. [3] by vertical 
dashed lines, which are obtained by requiring that the 



available volume fraction for arranging the conducting 
particles 

<Aavail. = (1 - <^2)('+"^/"^)' (1) 

coincides with the percolating volume fraction of voids, 
^$oid' which for the three dimensional penetrable sphere 
model used here is 0void — ^-^^ [^^J- 

As shown in Fig. [31(a), the slope of the initial decrease 
of if]\ is steeper for CF2/CF1 larger, and the position of the 
minimum gets shifted to higher values of 02- A similar 
effect is found by decreasing the penetrable shell thick- 
ness d for fixed (72 /ai. Fig. [31(b), leading to infer that 
for d/(j2 ^0 the minimum disappears and r]l decreases 
monotonously all the way up to 02- These features, and 
in particular the appearance of a minimum (i.e. opti- 
mal) value of the percolation threshold for finite pene- 
trable shells, represent our main finding and provide a 
previously unnoticed scenario for segregated percolation. 

Let us discuss now the physical origin of the non- 
monotonous behavior of the percolation threshold. The 
initial decrease of 77J can be fairly well reproduced by 
assuming that, for low values of 02, the volume frac- 
tion of the composite conducting particles (hard-core 
plus penetrable shell) is reduced by the volume occu- 
pied by insulating spheres. However, since the penetra- 
ble shells of the conducting particles may actually over- 
lap the insulating spheres, these latter may be treated as 
having effectively a smaller volume, ^eff ^ ^2, leading to 
0^(02) 02(O)(1 — 02'^eff/'^2)- Taking into account that 
insulating particles with (72 ^ a, where a is the mean 
distance between the closest surfaces of nearest neighbor 
conducting particles, should be ineffective in reducing 0^, 
we approximate t'eff by a sphere of diameter (72 — a. Fi- 
nally, by expanding 0i(02) in powers of /?i(02) — ^i(O), 
at the lowest order in 02 we find 

where 0i(O)' = lim02^o <^(/>i(02)/<^^i- As it is seen in 
Fig.[3l where Eq.([2j) (dotted lines) is plotted by using a = 
((71 +(i)/27^J(0)^/3 -(71 and 0i(O) as given in Ref.0, 
the low 02 behavior of rjl is rather well reproduced for 
all cases considered. 

By construction, the above argument neglects possible 
effects of 02 7^ on the connectivity number /c, i.e., the 
number of conducting particles directly connected to a 
given one. Actually, as it is shown in Fig. [T] where the 
color map defines k for each particle in the percolating 
cluster, the rather narrow /c-distribution for the homo- 
geneous case, which is peaked around the mean value 
{k) 2.25 [y, illl, changes drastically in the highly seg- 
regated regime of Fig.[Tl(b). Here, clusters of highly con- 
nected particles {k large) are bound together by "chains" 
of particles having low k values. Such distribution of k 
values is due to the fact that, in the vicinity of 02 7 
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FIG. 4: (Color online) Mean connectivity number (k) as a 
function of 02 for the same cases of Fig. [3ja). 



structure of the void space available for arranging the 
centers of the conducting particles is characterized by 
many narrow (quasi-one dimensional) necks connecting 
more extended void regions [20|. Percolation is possible 
only if such necks are populated by connected conduct- 
ing particles, and since for 02 ^2 necks become 
narrower, and so have less probability of being popu- 
lated, more particles are needed to ensure connectivity, 
thereby "overcrowding" the many void regions between 
the necks. The net effect of such mechanism, not cap- 
tured by Eq.([2]), is the enhancement of 77^ as 02 4^2- 
This is demonstrated in Fig. 2] where (/c), plotted for the 
same cases of Fig. [3l^a), displays a sudden enhancement 
(more marked for g^Icfx larger) at values of 02 corre- 
sponding to the points of upturn of v[{ of Fig. [31(a). The 
competition between the effect of reduced available vol- 
ume, which lowers r({ [Eq.([2])], and the enhanced connec- 
tivity at high segregation, which increases r]\^ is therefore 
at the origin of the minimum percolation threshold ob- 
served by us. 

Before concluding, let us discuss the possibility of ob- 
serving the features presented here in real segregated ma- 
terials. In conductor-insulator composites where trans- 
port is driven by tunneling, d represents the maximum 
tunneling distance between the conducting particles, so 
that d would be of the order of few nanometers. The same 
is expected for particles linked by organic molecules, as 
discussed in Ref.fil]]. For such values of the results 
of Fig. [3] would therefore apply to nano-composites with 
G\ d and G2 not exceeding a few tens of nanometers. 
Much larger values of d are however possible in some 
Ru02-glass composites, where a reactive layer of thick- 
ness 0.2-0.4 /im (or even more) surrounding the RUO2 
particles presents modified chemical and structural prop- 
erties [22| , most probably favoring hopping processes [23[ . 
In this case, the parameters used in our work would eas- 
ily account for composites with g\ in the range 50-500 
nm and G2 of few microns. 
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